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Abstract 

The effects of the CP violating 6 term in the QCD Lagrangian upon low energy 
hadronic phenomenology are reconsidered. Strong CP violating interactions among Gold- 
stone bosons and octet baryons are incorporated into an effective chiral Lagrangian frame- 
work. The 6 term's impact upon the decays rj nn and —>■ 77 is then investigated 
but found to be extremely small. A refined model independent estimate of nonanalytic 
contributions to the neutron electric dipole moment is also determined using velocity de- 
pendent Baryon Chiral Perturbation Theory. We obtain the approximate upper bound 
1^1 < 4.5 X 10-1°. 
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CP invariance is known to be a very good but inexact symmetry of the standard model. 
Violations of this discrete symmetry have been observed in neutral kaon decays and are 
hoped to be seen in bottom meson phenomena such as B — B mixing [0. CP violation 
in these weak processes may be attributed to a complex phase in the Kobayashi-Maskawa 
matrix. However, this phase is not the only source of CP violation within the minimal six 
quark standard model. In the strong interaction sector, instanton effects generate a GG 
term which preserves charge conjugation but breaks parity. This topological term enters 
into the QCD Lagrangian 



Cqco = -]Gii''G^,a+W<l-QMq + O^Gii''G^,a (1) 



with an undetermined coefficient 9 that represents a fundamental parameter of the standard 
model. In this letter, we investigate the effects of the 9 term upon low energy hadronic 
phenomenology. 

To begin, we restrict the QCD Lagrangian to three light quarks and take the mass 
matrix M to be real and diagonal without loss of generality: 




M = md . (2) 



It is convenient to rotate the 9 parameter away from the GG term in (|1|) and into the 
quark mass matrix which subsequently becomes complex : 

Cqcd -^G^^G^^a + qip(l -Wq - i9fnqj^q 

^ (3) 
= --GTG^.a + qJpqL + q^Tp qn -qAM + i9m)q^ -q^M- i9m)q^. 

If any of the current masses in (|^) equals zero, the reduced quark mass 

_ rriumdms 

m = 

ruumd + mdrus + msfriu 

vanishes and the mass matrix remains real. However in the real world, the current quark 
masses have the small but nonvanishing values 

(m„,md,m,) ~ (5, 9,181) MeV (4) 

which imply m ~ 3.2 MeV. So the 9 parameter disappears from the topological term's 
coefficient and reemerges in arg (det(M + i9m)) =9 + 0{9^). 
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At low energies, it is useful to replace Lagrangian (|^) which describes QCD in terms 
of fundamental quarks and gluons with an effective Lagrangian of mesons and baryons. In 
particular, the self interactions of Goldstone bosons associated with the chiral symmetry 
breakdown G = SU{3)r. x SU{3)n = SU{3)r.+ fl. may be analyzed in a nonlinear chiral 
Lagrangian framework 0. The Goldstone bosons appearing in the pion octet 
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(5) 



enter into the chiral theory through the combinations E = e^'^'^l ^ and ^ = e^"^! ^ where 
/ = 93 MeV represents the pion decay constant to lowest order. These exponentiated 
fields transform under G as 



(6a) 
(66) 



where L and R denote global elements of 5'C/(3)^ and 5'C/(3)h while local matrix \J{x) 
is implicitly defined by (^. Goldstone terms in the effective Lagrangian are constructed 
from E and ^ in a derivative expansion with respect to the chiral symmetry breaking scale 
^ 1 GeV. 

Incorporating baryons into this scheme seems problematic. The momentum expansion 
fails for baryons since their masses are not small compared to A^^.. However, a consistent 
approach for performing Baryon Chiral Perturbation Theory has recently been developed 
by Jenkins and Manohar f^^TlJ] using ideas and techniques familiar from the Heavy Quark 
Effective Theory |jl^,|13[. At low energies, a baryon may be regarded as an almost on-shell 
heavy particle that travels along a straight worldline. In this kinematic regime, a baryon 
has four-momentum p = Mj^v + k where the residual momentum k is small compared to 
its rest mass M^. As its four- velocity v is essentially unaffected by soft Goldstone boson 
absorption or emission, the baryon can be described by the velocity dependent field 

B^{x) = e*^«^-"S(x) 

which has the rest energy removed from its definition. The use of such fields allows one to 
formulate a well-behaved derivative expansion in terms of the small quantity k/h.^. This 
is the central idea behind velocity dependent Baryon Chiral Perturbation Theory. 
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The leading contributions to the effective Lagrangian 



V 

appear in d — 4 — e dimensions and at the renormahzation scale A as m 

= -^!-Tr{d^'E^ df,E) (7a) 
C'^\v) = iTvB^vVB, + I^M,7'^75{A^, B,} + FTrB„7'^75[A^, B,]. (76) 

Pions derivatively couple to the baryon octet through the Goldstone vector field 

that resides within the covariant derivative 

P'^B^ = a'^B, + [V^,B„]. 
In addition, they communicate via the axial current 

which explicitly appears in (ff^). We adopt the best fit values D = 0.61 and F = 0.40 
reported in refs. [1T0| , [TT[] for the symmetric and antisymmetric axial current couplings. We 



also include electromagnetic interactions into the effective Lagrangian by gauging a U{1)em 
subgroup of the global SU (3)^ x SU{3)u^ symmetry group and treating photons as external 



^ In their original formulation of Baryon Chiral Perturbation Theory, Jenkins and Manohar 
introduced spin operators 5^ that act on the velocity dependent baryon fields and incorporated 
them into the chiral Lagrangian. For the applications that we will consider, these spin operators 
offer no real advantage over conventional gamma matrices. We therefore express the derivative 
interactions of Goldstone bosons with baryons in our Lagrangian in terms of the latter rather 
than former objects. 

^ The baryon decuplet may also be readily included into the low energy Lagrangian p|-pT|. 
However, these fields do not influence CP violating baryon octet phenomena such as the neutron 
electric dipole moment at leading nontrivial order. We consequently neglect them here. 



3 



fields. The partial derivatives acting on hadron fields in (|^ are then promoted to covariant 
derivatives with respect to electromagnetism 

d^B, d^By - iA'/^eA^[Q, B„] 

where 



Q 



denotes the quark charge matrix. 

The symmetry breaking efi^ects of the quark mass matrix = M+i9fn in the underly- 
ing QCD Lagrangian can be systematically incorporated into the low energy theory via the 
spurion procedure. We first regard as a fictitious field and assign it the transformation 
rules 

M^LMR^ 

which preserve chiral symmetry and discrete CP in eqn. (|^) . We then construct G and CP 
invariant terms in the chiral Lagrangian from the hadron fields and matrix Ai. Finally, 
we set the spurion to its true constant value and thereby generate interactions that violate 
chiral symmetry, fiavor and CP in the effective theory. 

Working to linear order in Ai, we add the following terms into the low energy La- 
grangian: i 

^(M) ^ ^L!Z!^Tr(EA^t + ME^) (8a) 
C'^^\v) = 6^TrB„{^t_^^t + ^_^t^^ + 6^TrB„[^t_^^t + ^_^t^^ b„] 

+ (7Tr(^(E-l)A^^ + A^(E^-l))TrB„B„. (86) 

Decomposing the mass matrix into its real and imaginary parts, we obtain CP conserving 
terms proportional to M 

^iM) ^ ^L!Z!^TrM(E + Et) (9a) 
J^'s'^v) = 6^TrB4^tM^t + ^m^, BJ + 6^TrB„[^tM^t + B„] 

+ (TTrM(E + E^ -2)TrB„B„ (96) 



^ A common baryon mass shift has been extracted from the a term in (86) and implicitly 
absorbed into the average octet mass parameter Mb- 
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and CP violating interactions linear in Om 

= _^__^^6'mTr(E - E^) (10a) 

C^^\v) = -iem\bnTrBy{E - E^ B4 + 6^TrB^[E - E^ B^] + aTr{E - E^)TrB^B„ . 

(106) 

The leading terms in ( |9a, l\ ) produce Goldstone boson masses and split the baryon octet 
multiplet. The parameters and bp as well as the mean baryon mass Mg can therefore 
be fixed by a fit to the hadron mass spectrum using the assumed current quark masses in 
(p. We find the approximate tree level values 

a = 1302 MeV = 1197 MeV 

(11) 

bn = 0.18 bp = -0.54. 

Having set up the necessary formalism, we can now investigate specific examples of 
strong CP violating phenomena. We first consider Goldstone boson processes. Expanding 
pion Lagrangian ( |10a|) , we isolate the self interaction terms 

J^'^' = -^^/x^rfrrr(7r3) +0(7r5) 



/ 



IJ,9m 



2v3 , _ vS n r, v3 I v3 ^^n-TT^ v3 o 
?77r+7r H t^tt^tt" r]K^K 'qK^K r]-^ 



3 ' 3 ' 3 ' 9 



+ 0(7r5). 
(12) 

Prom these trilinear vertices, one can easily compute the tree level rate for the CP violating 
decay ij tttt P,!!^: 



1271 V / / 

Comparison with the current experimental limit r(r^ '^~^'^~)exp < 1-79 x 10~^ MeV 
implies that \9\ < 4.6 x 10~^. However, much more stringent bounds on 9 are available from 
electric dipole measurements. So it is of greater use to insert \9\ < 4.5 x 10~^°, the con- 
straint which comes from the neutron electric dipole analysis that we will shortly present, 
into eqn. (^. We then obtain the prediction T{r] n~^'n:~)th < 1-7 x 10~^° MeV. This 
theoretical bound is unfortunately 14 orders of magnitude smaller than the experimental 
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upper limit. One may thus safely assume that the 9 parameter will never be extracted 
from this Goldstone decay mode. 

The CP violating interactions in (|12D also contribute to loop processes. The simplest 
one-loop diagrams involving the Goldstone vertices are the tt^ and ry tadpoles displayed in 
fig. |l|. The graphs' 0{M.) quadratic divergences automatically sum to zero as a result of 
the vanishing tadpole renormalization condition 

(0|z^mg7^g|7r") = 

which we implicitly invoked when rotating 9 into the quark mass matrix . Terms that 
persist at higher orders in Ai along with tadpole interactions in the effective Lagrangian 
starting at 0(A^^) may always be eliminated via an appropriate modification of the 9 
rotation. 

While neutral Goldstone bosons do not disappear directly into the vacuum, they can 
decay into two photons. As illustrated in fig. H, the process tt^ 77 may proceed through 
kaon triangle and "seagull" graphs that violate CP. A straightforward calculation yields 
the finite and gauge invariant on-shell amplitude 



2fx9fne 



-—- l^ltan , = -1 [g^^-2^f)e^{p)e [q). 



(14) 



From the form of this expression, we see that the graphs in fig. ^ match onto an infinite 
string of nonrenormalizable operators that reduce via the equations of motion to ti^F^^ F^jj. 
Since this operator does not respect CP, it cannot interfere with the one pion-two photon 



vertex in the Wess-Zumino action 17,18 



W{T.,An = I d'^{-^^'F^''F,. + ---}- (15) 
The CP violating contribution to the decay rate 

r(.« r,)^ = ±. (^)\^Vm-^ (tan- ^ , ^ Y - l] ' (16) 

647r \ n J \ f J I mi \ Am^r. - / J 



is therefore suppressed by two powers of 9 compared to its lowest order CP conserving 
counterpart 
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Their ratio is consequently miniscule: 

r(7r0 77)^2 



= 2.1 X 10"^^^ < 4 X 10"^^! 



The preceding Goldstone boson examples clearly demonstrate that one cannot hope 
to observe strong CP violating effects unless they occur at linear order in 9. Fortunately, 
there is a well-known phenomenon in the baryon sector which meets this criterion: the 
neutron electric dipole moment (NEDM). Previous attempts to determine the magnitude 
of this important observable have typically relied upon specific model calculations. We will 
reconsider this problem in the context of velocity dependent Baryon Chiral Perturbation 
Theory. While our investigation is clearly similar in spirit to the well-known current 
algebra NEDM analysis of Crewther, Di Vecchia, Veneziano and Witten 0, we believe 
that the effective field theory approach is more transparent and systematic. Furthermore, 
heavy hadron techniques are especially well suited for studying the electric dipole moment 
(ijv since it is a static property. Recall that is defined in terms of the form factor F3 
appearing in the interaction Lagrangian 

^' = ^^^N{p')a^'^^'N{p)F,^ (18) 

evaluated at zero momentum transfer [q = p' — p = 0) : 

Fs{0) 



dp 



The neutron's four-velocity is therefore conserved, and the assumptions underlying the 
static hadron picture are genuinely satisfied. So we view the NEDM problem as a partic- 
ularly nice application of Baryon Chiral Perturbation Theory. 

We first enumerate the leading terms which may appear in the low energy Lagrangian 
and directly contribute to dpf. Group theory counting indicates that there are ten indepen- 
dent terms which one can form from the octet fields By and By, the photon combination 
eQa^^'-y^F^jj, and the parity-odd mass term ^^A^^^ — ^A^^^. A convenient basis for these 
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dimension-six operators is given below: 



Oi = 




02 = 


-^Tr(B„a'^-75F^4Q,B„])lV(A^Et - EA^t) 


Os = 


-^Tr(B„a'^-7'i^M.B„)Tr(Q(e^A<^^ - ^A^^O) 
^x 


04 = 


-^IV(B„a'^-7^F^4et_^^t _ ^A^t^))^(QB„) 
^x 


05 = 


-^Tr(B„a'^-75F^.Q)Tr((^t^^t _ ^M^^)B,) 




-^IV(B„a'^Vi^^.[(?W -^-M^?)Q,B.]) 
^^x 


O7 = 


4rTr fB„(T''^7^F„^ [Q, B„l {(^ Mi^ - f Al V)) 

X 


Os = 


-^IY(B„a'^-7^F^4Q, (e^At^t _ ^AltOB„]) 

X 


09 = 


-^(Tr(B„a^Vi^M-B„[^t^^t _^_Mt^,Q]) 

X 




+ Tr(B„a^-75F^4^t^^t _ ^^t^^ qb„])) 



Of these ten operators, only Oi affects the electric dipole moment of neutral baryons. Its 
impact upon is unknown however since Oi appears in the effective Lagrangian with an 
a priori undetermined coefficient. In principle, the coefficients of all the operators in (|T9| ) 
could be fixed from Goldstone boson-photoproduction data. But in practice, the challenge 
of performing such a fit appears formidable. 

As first noted in ref. 0, the most important contributions to the neutron's dipole 
moment actually do not come from tree level composite operators but rather from one- 
loop Goldstone boson diagrams. Such graphs generate infrared log renormalizations of Oi 
which diverge in the chiral limit. Their nonanalytic dependence on M. is exactly calculable 
unlike the analytic dependence of Oi's coefficient. The nonanalytic terms consequently 
provide a rough but useful order-of-magnitude estimate for d^- 

There are a number of CP violating diagrams that enter into the neutron-neutron- 
photon IPI Green's function r*^^^'*'^ at one-loop order. But only those listed in fig. ^ yield 
chiral log corrections to the NEDM. In addition, these graphs contain fractional power 
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terms starting at 0{Ai^^'^) which result from the mass sphtting between the external neu- 
tron and intermediate baryon. Summing the diagrams with internal tt" and P propagators 
as well as and E~ lines, we obtain 



167r2/2 



1 

log — 



7r(Mp - Mj, 



log 



A' 



m 



K 



~{D-F){hD-hF) 
The Dirac algebra identity 



ml - {Mp-M^Y^ 
7r(ME - Mn) 



K 



(Me - Mat) 2 



(20) 



transforms the velocity dependent expression into a manifest dipole operator. Numerically 
evaluating eqn. (|0D at the scale A = = 1 GeV, we find for its coefficient 



(2.76 - 0.08) X 10-^^^e-cm iiV^(A;')a'^V(^' - k)^N^{k) + 



where the pion and kaon contributions are separately displayed. The large disparity be- 
tween these terms stems in part from their SU{3) couplings and infrared logarithms. But 
more importantly, the discrepancy results from the near cancellation between the logA4 
and terms in the strange virtual hadron graphs. The total nonanalytic contribution 

to the neutron electric dipole moment is thus given by 

= (2.68 X 10"^*^^)e-cm. 

Comparing with the current experimental upper limit {d^l < 1.2 x 10~25g.c]2i at 95% CL 
111], we deduce \e\ < 4.5 x 10"^°. 

In conclusion, the results from our chiral Lagrangian investigation of the 9 term's effect 



upon low energy hadron phenomenology are in basic accord with earlier findings Its 
virtue is therefore not novelty but rather simplicity. Of course, as in any effective field 
theory analysis, the results are model independent and can be systematically improved 
by retaining higher order terms in the derivative expansion. In particular, nonanalytic 
0{1/Mg) corrections to the NEDM could be determined. However given the uncertainty 
in the analytic contributions to |(ijv|, such further refinement of our estimate for 9 is not 
of much practical importance. 
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Added note: After completion of this work, we learned that results similar to those pre- 



sented here have been reported in ref. [21]. See also ref. [22 
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Figure Captions 



Fig. 1. One- loop tt and ry tadpole diagrams. Large dots denote CP violating vertices. 

The cross in the tv^tv^t] graph represents the OirUu — TUd) off-diagonal mixing 
term in the neutral Goldstone boson propagator. 

Fig. 2. CP violating graphs that contribute to 7r° — > 77. A virtual charged kaon runs 
around the loop. 

Fig. 3. One-loop diagrams that contribute to the neutron electric dipole moment. Iden- 
tical graphs with internal 7r~ and Py propagators replaced by K'^ and E~ lines 
are not pictured. 
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